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Abstract 

The recent measurements of azimuthal single spin asymmetries by the HERMES collab- 
oration has opened up new possibility to measure chiral-odd distribution functions through 
semi-inclusive deep-inelastic scatterings. Here, predictions are given for the twist- 2 and twist- 
3 chiral-odd distribution functions of the nucleon within the framework of the chiral quark 
soliton model, with full inclusion of the vacuum polarization effects as well as the subleading 
1/N C corrections. The importance of the vacuum polarization effects is demonstrated by 
showing that the so-called Softer inequality holds not only for the quark distributions but 
also for the antiquark ones. 



It is widely believed that the transversity distribution function hi(x) provides us with 
valuable information for our thorough understanding of internal nucleon spin structure 
Because of its chiral-odd nature, however, it does not appear in the standard deep-inelastic 
inclusive cross sections (at least in the chiral limit). The only practical way to determine it has 
been believed to use Drell-Yang processes. Recently, the HERMES collaboration proved for the 
first time that it is practically feasible to measure chiral-odd distribution functions by making 
use of the so-called Collins mechanism in the semi- inclusive electro-pion productions ||. 
The first HERMES measurement of the azimuthal single target-spin asymmetries for charged 
pions has been done by using longitudinally polarized proton target. From the theoretical 
viewpoint, much simpler would be the measurement of the similar asymmetry for the target 
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polarized transversely to the incident electron beam, since it enables us to determine the 
transversity distribution hi (x) more directly and efficiently Such measurement will in fact be 
performed soon. 

In view of this very exciting experimental situation, it is highly desirable to give some 
reliable predictions for the chiral-odd distribution functions of the nucleon. Naturally, complete 
understanding of the parton distribution functions (PDF) need to solve nonperturbative QCD 
dynamics. Unfortunately, the direct evaluation of the PDF in the lattice QCD simulation is 
not possible yet. At the present moment, we must therefore rely upon some effective theory 
of QCD. We would claim that, among many of them, most promising would be the chiral 
quark soliton model (CQSM) 0-0. In fact, a prominent feature of the CQSM is that it 
can simultaneously explain two biggest findings in the recent experimental studies of high- 



energy deep-inelastic scattering observables, i.e. the EMC measurement in 1988 [|Kj, which 
dictates unexpectedly small quark spin fraction of the nucleon, and the NMC measurement in 
1991, which has revealed the excess of d sea over the u sea in the proton JTTJ] . Furthermore, 
its unique prediction, i.e. the flavor asymmetry of the spin-dependent sea-quark distribution 
Au(x) — Ad(x) > JT2] , , seems to win some semi-phenomenological and/or semi-theoretical 
support lfL4|1 ) although not entirely definite yet. 

What is the chiral quark soliton model, then? For large enough N c , a nucleon in this 
model is thought to be a composite of N c valence quarks and infinitely many Dirac sea quarks 
bound by the self-consistent pion field of hedgehog shape 0. After canonically quantizing 
the spontaneous rotational motion of the symmetry breaking mean field configuration, we can 
perform nonperturbative evaluation of any nucleon observables with full inclusion of valence 
and deformed Dirac-sea quarks . This incomparable feature of the model enables us to make 
a reasonable estimation not only of quark distributions but also of antiquark ones, as we shall 
demonstrate later. Finally, but most importantly, only 1 parameter of the model (dynamically 
generated quark mass M) was already fixed by low energy phenomenology, which means that 
we can give parameter-free predictions for parton distribution functions at the low renormal- 
ization scale. There already exist some investigations of the chiral-odd distribution functions 
based on the CQSM (or the Nambu-Jona-Lasinio soliton model). The first investigation of the 



transversity distribution hiix) in the CQSM was done by Pobylitsa and Polyakov [15|. Al 



though pioneering, their investigation should be taken as qualitative in several respects. First, 
the effect of vacuum polarization is entirely neglected. Second, only the isovector combination 
of hi{x), i.e. hi(x) — hf(x), was estimated, and the isocalar piece h^(x) + hf(x), which arises 
only at the sub leading order of 1/N C expansion, was left over. On the other hand, Gamberg 
et al. investigated both of h^ d (x) and h\ (x), but also within the "valence-quark-only" ap- 
proximation ||16|| . We shall later show that this approximation does not make full use of the 
potential power of the model and cannot be justified. The purpose of the present investigation 



is to give more complete predictions for the chiral-odd distribution functions with full inclusion 
of the vacuum polarization effects as well as the sub leading 1/N C corrections. 

We start with the standard definition of the chiral-odd spin-dependent PDF H : 

h[ I=0/I=1) (x) = hl{x)±h d Ax) = -L f° dz e ixMNZ0 

4 71 J-oo 

x (iV(P = 0)5 ± |^(o)(l + 7 V) 7±75 I M i,(z)\N(P = 0)S 1 _)\ za= . Z(hZ±=0 , (1) 
h ( l =0/I=1) (x) = h u L {x)±h d L (x) = — [°° dz e txMNZ0 

ATT J-oo 

x (N(P = 0)S Z | V f (0) 7 3 75 | * | I N(P = 0)S Z ) \ Z3 =- Zo , z±=0 . (2) 

What we must evaluate here is the nucleon matrix elements of quark bilinear operators with 
light-cone separation. As fully discussed in f|,|§, on the basis of the path integral formulation 
of the CQSM, such nonlocality effects in time and spatial coordinates can be treated in a 
satisfactory way. Omitting the detail, we just recall here the fact that, within the framework of 
the CQSM, the isoscalar and isovector distributions have totally dissimilar theoretical structure 
since they have different dependence on the collective angular velocity Q : 

h u 1/L (x) + h d 1/L (x) ~ N c O(^) ~ 0(N° C ) , (3) 

K /L (x) - h d 1/L (x) ~ N e [0(Nl) + Oin 1 )] ~ 0{Nl) + 0(N° C ), (4) 

where use has been made of the fact that Q scales as 1/N C . A noteworthy feature here is the 
existence of the subleading 1/N C correction to the isovector combination of the distributions. Its 
importance is already known from a similar analysis for the isovector longitudinally polarized 
distribution M or its first moment, i.e. the isovector axial coupling constant of the nucleon 

n • 

We summarize in Fig.l the theoretical predictions for the chiral-odd distribution functions 
hi(x) and h^x). Fig. 1(a) and Fig. 1(b) respectively stand for the isoscalar and isovector 
combinations of /i"(x) and hf(x), while (c) and (d) represent the similar combinations of h^(x) 
and h d L (x). Here, the distribution functions with negative x should be interpreted as aitiquark 
ones according to the rule : 

h u 1/L (-x) ± h d /L (-x) = - [h* /L (x) ± h d /L (x)] (0<x< 1). (5) 

In all the figures, the long-dashed curves peaked around x ~ 1/3 represent the contributions of 
N c valence quarks, while the dash-dotted curves are those of Dirac-sea quarks in the hedgehog 
mean field. The sums of these two contributions are denoted by solid curves. Concerning the 
transversity distributions h\(x), one sees that the effects of Dirac-sea quarks are not very large. 
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One also notices that both of h\{x) + hf(x) and h^(x) — hf(x) have fairly small support in 
the negative x region, which means that the transversity distributions do not have significant 
antiquark components, in contrast to the unpolarized as well as the longitudinally polarized 
distributions. Turning to the twist-3 distribution Iil(x), one finds that the effect of Dirac-sea 
quarks are sizably large. Nevertheless, an interesting feature is that a considerable cancellations 
occurs between the contributions of valence quarks and of Dirac-sea quarks in the negative x 
region. As a consequence, the antiquark components are not so significant also in the case of 
h L (x). 
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Figure 1: Theoretical predictions for the chiral-odd distributions hi(x) and hi(x). (a) and 
(b) respectively represent the isoscalar and isovector parts of h\(x), while (c) and (d) are the 
isoscalar and isovector parts of hi(x). In all the figures, the long-dashed and dash-dotted 
curves denote the contributions of N c valence quarks and of the Dirac-sea quarks, whereas the 
solid curves are the sum of these two contributions. 

The transversity distribution h±(x) is of special importance, since it forms the set of twist-2 
distribution functions, together with the unpolarized distributions j\ [x) and the longitudinally 
polarized ones g\{x). It is known that these distribution functions must satisfy the so-called 
Soffer inequality |18| expressed in the following form : 

\h\{x)\ < \ {±ft{x) + 9l(?)) (x>0,x<0). (6) 



Here, the plus sign of ff(x) corresponds to the region x > 0, while the minus sign to 
x < 0. (Note that the inequality with negative argument gives a relation for the antiquark 
distributions.) Now the question is whether the predictions of the CQSM fulfill this general 
inequality or not. Fig. 2 shows that, if one includes the vacuum polarization effects properly, the 
Soffer inequality is well satisfied for both of the w-quark and the <i-quark. On the other hand, 
if one neglects the Dirac-sea contributions, it is badly broken for the antiquark distributions. 
An important lesson learned from this observation is that the field theoretical nature of the 
model, i.e. the proper inclusion of the vacuum polarization effects, plays essential roles in giving 
reasonable predictions for antiquark distributions. Another lesson is that the frequently-used 
saturation Ansatz of this inequality for estimating h\{x) is not justified. This seems reasonable 
to us, since the magnitude of ff(x) is much larger than g\(x) and since h\(x) is closer to g\(x) 
rather than ff(x). 
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(b) with Dirac-sea quarks 
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(c) valence-only approx. 
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(d) valence-only approx. 
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Figure 2: Theoretical check of the Soffer inequality, (a) and (b) are the results obtained with 
full inclusion of the vacuum polarization effects, while (c) and (d) are obtained with "valence- 
quark-only" approximation. In all the figures, the solid and dashed curves represent the r.h.s 
and the l.h.s. of the inequality (6), respectively. 



As already touched upon, HERMES group recently carried out a very interesting measure- 
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ments on single target-spin asymmetries in semi-inclusive pion electroproduction ||. What 
they have measured are the azimuthal asymmetries, i.e. the asymmetries of the cross sections 
depending on the azimuthal angle <fi with respect to the lepton scattering plane : 



.4 



w 

UL 



J d(j)dy W(4>) (da + /S^ dxdydcf) — da /S H dxdydcf)) 



(7) 



| J d(f>dy (da+fStf dxdydcf) + da~/S H dxdydcf)) 
with W(cf>) = sin0 or sin 20 and with S H being the nucleon polarization. The theoretical 



analyses of these azimuthal asymmetries were already performed by several authors |19[-p3 
Here we use the expression given in EI : 
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(12) 



withx = Q 2 /2(P-q),y = (P-q)/(P-k 1 ) and z = (P-P h )/(P-q) in the notation of E3|. One re- 



alizes that the asymmetries depend on 4 distribution functions fi(x), hi(x), h L (x), hf^\x) and 
3 fragmentation functions Di(z), H± (z), H(z). Here D\(z) is the familiar spin- independent 
fragmentation function, while Hi^ l \z) is (/cy/2M|)-weighted fragmentation function defined 
by 

\k T ? 



Ht (1)a (z) 



= d 2 k' 



T 



2 Ml 



Ht a (z,k' T ) 



(13) 



with k' T = —zkT, where H^~ q (z, k' T ) is a T(time-reversal)-odd leading twist fragmentation 
function, giving the probability of a spinless or unpolarized hadron to be created from a 
transversely polarized scattered quark. The remaining subleading T-odd fragmentation func- 



tion H a (z) is known to be constrained by the relation H a (z) 



.(l)a, 



due to the 



Lorentz invariance [^0j. To estimate these T-odd fragmentation functions, we follow [22|-[23] 
and use the Collins Ansatz 101 



A c {z,k T ) = 



\kx\ H^~(z,k 



M h 



Di(z, k 2 T ) 



M c \k T \ 
M 2 +k, 



2 ' 

T 



(14) 



with the parameters Mq = 2 m n and rj = 1.0. 

Coming back to the distribution functions, we need 4 functions fi(x), hi(x), and 
^il\ x )- I n the following, simply by using the GRV parameterization of the unpolarized 
distribution fi(x) p4"l , we concentrate on the chiral-odd distributions. Here, h\(x) is the 
familiar transversity distribution function. The function hi(x) consists of the twist- 2 part, 
which is expressed by hi(x), and the interaction dependent part or the genuine twist-3 part, 
as 

h a L (x) = 2x f 1 dy^f- + h a L (x). (15) 
Jx y 2 

On the other hand, the third function (x) is expressed with use of hf(x) and h a L {x) as 

htP\x) = f dy [hl(y) - ht(y)]. (16) 

In the absence of powerful theories, the following two approximations have been frequently 
used for estimating these functions. The first is the twist-2 or the Wandzura-Wilczek type 
approximation for h a L {x), i.e. the approximation to set h a L {x) ~ 0. The second approximation 
is to assume approximate equality of h\{x) and h1(x), which leads to hf^ a (x) ~ 0. This 
second approximation is advocated by HERMES group, by the reason that it seems consistent 



with small sin 20 asymmetry obtained by their experiment [22]- [23]. 

Since the CQSM can give some reasonable predictions for both of hi(x) and h^x), the 
other distributions can also be determined by using (16). Fig. 3 summarizes the theoretical 
predictions of the CQSM for these chiral-odd distributions with different flavors. The distri- 
butions shown here corresponds to the energy scale of Q 2 = 2.5 GeV 2 , an average value for the 
HERMES kinematical region. The scale dependencies of the distribution functions are taken 
into account in the following way. First, to include the scale dependencies of hi(x) and the 
twist-2 part of hi(x), we use the Fortran program of leading-order DGLAP equation |2~5| . On 
the other hand, the scale dependence of the twist-3 part of hi(x), i.e. hi(x), is taken into 
account by solving the leading-order DGLAP type equation obtained in the large N c limit 
|2(| , [^7] . The starting energy of these evolutions is taken to be Q 2 = 0.23 GeV 2 [|28|| . From 
Fig. 3 (a) and 3(b), one reconfirms sizable differences between the two distributions hi(x) and 
^l(x), especially at lower values of x. A general tendency is that Iil{x) is concentrated in the 
lower x region as compared with h\{x). We recall that this characteristic is also observed in 
the MIT bag model predictions for the same distributions. In fact, the chiral-odd distributions 
hi(x) and Hl{x) in both models are not extremely different as far as the dominant u-quark 
distribution in the proton is concerned. (This is not necessarily true for the corresponding 
(i-quark distributions as well as the u and d ones. We also recall that the situation is quite 
different for the other two twist-2 distributions, i.e. the longitudinally polarized one and the 
unpolarized one.) The sizable difference between h\(x) and hi{x) indicates that the Ansatz 
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hi(x) ~ Hl{x) adopted by the HERMES group may not be necessarily justified. To verify it, 
we plot in Fig. 3(d) (x) obtained from (16). One sees that (x) for u-quark has in fact 
nonnegligible magnitude, although the distributions for d, u and d are pretty small. We also 
show in Fig.3(c) the twist-3 part of h^ix) obtained with (15). They are generally small but 
certainly nonnegligible especially in the small x region. 
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Figure 3: Theoretical predictions for the chiral-odd distribution functions h 1 (x),hL(x),hi J (x) 
and hi^\x) for each flavor at the energy scale Q 2 = 2.5 GeV 2 . 



Fig. 4 shows a preliminary comparison between the predictions of the CQSM and the HER- 
MES experiments for the sin0 and sin 20 asymmetries. The theoretical curves have been 
calculated as in [22|,|23| by integrating over y and z in the HERMES kinematical range with 
the parameter a = 0.44 GeV and b = 0.36 GeV by assuming Gaussian distribution for the 
distribution functions and the fragmentation functions. The main difference with the similar 
analysis done in the same model is that we have included all the subleading contributions 
into the distribution functions (no twist-2 approximation) as well as into the fragmentation 
functions. In fact, it turns out that the subleading term containing H a (z) in (11) gives sizable 
contributions to the sin asymmetry. One can say that the theory reproduces the qualitative 
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Figure 4: Theoretical predictions of the CQSM for the sin0 and sin 20 asymmetries in the 
semi-inclusive electro-pion productions in comparison with the HERMES experiment. 

features of the data, although the uncertainties of the present experimental data are still too 
large to draw any decisive conclusion. As pointed out before, the HERMES group advocates 
that the observed small sin 20 asymmetry is consistent with the Ansatz h\{x) ~ Hl{x), or 
equivalently hf^(x) ~ 0. (Note that the formula for sin 20 asymmetry is proportional to this 
distribution function.) On the other hand, the CQSM as well as the MIT bag model indicates 
that these two functions are rather different especially in the small x domain. Unfortunately, 
it seems difficult to draw a definite conclusion only from the present HERMES data for the 
sin 20 asymmetry. We certainly need more accurate experimental data. 

From the theoretical viewpoint, the asymmetry obtained with a target polarized trans- 



versely to the direction of incident electron beam is much simpler |[20|| . Such asymmetry is 
proportional to 

Aot' ~ ^r 1 \St\ (1-y) £ 4xht(x) Ht Wa {z). (17) 

V a 

By measuring it, we can then get more direct information on the transversity distribution 
hi(x). Under the dominant-flavor-only approximation for the fragmentation functions, the 
semi-inclusive tt + and 7r _ productions respectively measure the following combinations of the 
transversity distributions : 

7T+ : t h\{x) + 1 hf(x), (18) 
7T- : I h*(x) + I h*(x). (19) 
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Figure 5: Theoretical predictions for the combinations of the transversity distributions, 
^hi(x) + ^hf(x) and \hf(x) + ^h^(x) at Q 2 = 2.5 GeV 2 , which are to be measured in semi- 
inclusive 7r + and 7r~ productions with a target proton polarized transversely to the direction 
of incident electron beam. 

The solid curves in Fig. 5 stand for the full predictions of the CQSM for these distributions 
at Q 2 = 2.5 GeV 2 , while the dashed curves here are obtained by dropping the antiquark 
components in these combinations. One sees that the first combination is quite insensitive to 
the presence of the antiquark component, while the second one is sensitive to the antiquark 
component at least in the small x region. This indicates that ir~ semi-inclusive production 
may be a useful tool to probe the role of antiquark in the transversity distributions. 

In summary, we have given theoretical predictions for the chiral-odd distribution functions 
h\{x) and hi{x) within the framework of the CQSM, with full inclusion of the vacuum po- 
larization effects as well as the sub leading 1/N C corrections. The importance of the vacuum 
polarization effects has been demonstrated by showing that the so-called Soffer inequality holds 
not only for the quark distributions but also for the antiquark ones. The theoretical predic- 
tions of the model for the azimuthal single spin asymmetries in the semi-inclusive electro-pion 
productions are shown to be qualitatively consistent with the corresponding HERMES data, 
although we certainly need more accurate experimental data as well as more precise knowledge 
about the T-odd fragmentation functions before drawing any decisive conclusion. We also find 
that the chiral-odd distribution functions have rather different magnitudes and a;- dependencies 
as compared with the longitudinally polarized ones investigated in previous works. Especially 
interesting is that the transversity distribution have very small antiquark components. We hope 
that these unique predictions of the CQSM will be tested through near future semi-inclusive 
meaurements, which enables the flavor as well as the valence plus sea-quark decompositions. 
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